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ABSTRACT 

The absolute magnitude of the brightest galaxy of clusters varies remarkably 
little and is nearly independent of all other physical properties of the cluster 
as, e.g., its spatial extension or its richness. The question arises whether the 
observed small scatter is compatible with the assumption of dynamical evolution 
of the cluster. This is investigated with the help of statistical analysis of the 
results of cluster simulation. The underlying interaction process is merging (and 
also destruction) of smaller galaxies forming the giant galaxy. The cluster itself 
is supposed to be in virial equilibrium. 

We find that the evolutionary importance of merger processes grows with 
decreasing scale. Rich clusters as well as their brightest members evolve merely 
slowly whereas compact groups as well as their brightest members evolve more 
rapidly and more violently. We also find that the number of merger processes 
leading to the growth of the brightest cluster galaxy (BCG) is small enough 
to keep the BCG mass dispersion below the measured value. Our simulation 
substantiate that just the combination of the initial distribution function and the 
following merging to form the BCG can explain the remarkably small variance 
of mean BCG masses between clusters of different size and different number of 
galaxies. 

Subject headings: Galaxies: clusters: general - Galaxies: evolution - Galaxies: 
interactions - Galaxies: mass function 
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Introduction 



The mean absolute magnit ude of brightest cluster gal a xies is {Myi) = —22. 7 mag with a 
dispersion of (t(Mvi) = 0.4mag (IHoessel. Gunn. &: ThuanI (llQSOl )). The di stribution of BCG 



absolu te magnitudes is well represented by a Gaussian distribution (e.g. 
(119951 ) ) , and the relative variation of luminosities yields 

a{L) 



Postman &: Lauei 



L 



0.5. 



Batcheldor et al.l (120071 ) show that the luminosity dispersion in the near infrared is even 
less. Due to th is small i ntrinsi c dispersion BCGs are ide al tools for cosmological distance 
determination (jSandagd (119721 ): iPostman fc Laued (119951 )). Furthermore their theoretical 
investigation offers the possibility to understand the properties and evolution of the cluster 
itself and its members. 



Sandagd (119721 ) points out that there is a variety of possible accounts for the observed 



scatter in BCG apparent magnitudes, but all contributions but intrinsic dispersion (i.e. 
scatter in absolute magnitude) seem to remain small, in particular is there no systematic 
trend with redshift. At least the small total dispersion shows that any single dispersion 
contribution remains small. 

There are many kinds of suggestions trying to explain the observed small scatter in 
magnitudes of brightest cluster galaxies. One possibility is the existence of a common lu- 
minosity function all galax ies of a cluster (including the brightest one) are drawn from 
( ISchechter fc Peebles! ( 119761 )). either with or without further development. The alternative 
consists of a unique creation process, a "standard mold", or a not yet discovered phys- 
ical effect that delimits the luminosities of galaxies and cuts off the bright end of the 
luminosity function. This hypothesis is underpinned by the astonishing fact that bright- 
est galaxies of rich clusters are just as luminous as the brightest observed field galaxies 
( iHumason. Mayall. fc Sandagd (119561 )). 



How special are BCGs? It has been argued that surface brightness profiles of BCGs are 
well fit by the same Sersic law that describes less-luminous spheroids and o bey the same rela- 



tions b etween fitting parame ters that characte r ise E/ SO galaxies generally (IBatcheldor et al. 



( 120071 )). On the other hand, Iv.d. Linden et al.l (120071 ) show that BCGs contain a larger frac- 
tion of dark matter as compared to non-BCGs of the same stellar mass, influenced by the 
cluster environment. 



X-ray observations of the intracluster medium (e.g., iPeres et al.l (Il998l )). in particular 



of the cluster core, has given rise to the assumption of c ooling flows depos iting up to 100 Mg 



yr ^, building up the BCG over many Gyrs. However, iMotl et al.l (120041 ) support the alter- 
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native that cores of cool gas are being built from the accretion of discrete stable subclusters. 
At the same time they realise that the common presence of substructure in galaxy clusters 
argues for accretion or m erger e vents that occured sufficiently recently not to be erased by 
relaxation of the cluster ( iBirdI (Il994l )). 



Galactic cannibalism, i.e. the accretion of the existing galaxy population through dy- 
namical friction and tidal strip ping predominantly in the e volved cluster, has been pro- 
posed by severa l auth o rs (e.g., lOstriker &: HausnianI (19771 )). Further in vestigat i ons b y 



MerrittdMerrittI fll983r i. iMerrittI (ll98J) 



Merritt 



(119851 )) and others (e.g., iLaueii (119851 ) 



Dubinski Jl998h ) showed, however, that galaxies are moving too fast and are significantly 
tidally truncated thereby suppressing dynamical friction. They concluded that the dynami- 
cal friction timescales are generally too long such that galaxies could therefore only accrete 
a small fraction of their total luminosity within a Hubble time. 

Galactic mergers remain as the most likely dynamical effect that dominates BCG evo- 
lution. But given a vivid merger history, it is problematical how a uniform magnitude of the 
brightest cluster galaxy can get over frequent encounters while being exposed to mergers with 
other galaxies, or deletion. Tightly related to this puzzle, the question comes up whether 
the common luminosity function of all up to date analysed clusters is a direct consequence 
of the initial mass fiuctuations of the early universe - or otherwise, if the mass spectrum at 
that time was an entirely other one, and evolution processes provide for the universality of 
the luminosity function. 



However, the scenario seems less agitated. iHoessell (119801 ) estimates, correlating spatial 
extension and luminosity of BCGs, that mergers of galaxies of comparable size lead to the 
formation of the first-ranked member. As he divides the collapse time of a cluster by the time 
that pass between two close encounters of two relatively big galaxies, he obtains about four 
merger events forming the BCG. This result coincides with the number of observed cores in 
cD-galaxies. In the same way, the luminosity function chan ged only little, because merger 
proces ses operate effectively only on large massive galaxies (ITrevese. Cirimele. &: Appodia 
Jl996h l 



Numerical simul ations of large groups of 50-100 spheri c al galaxies in an equilibrium clus- 



ter ha ve been done by lFunato. Makino. fc Ebisuzakil (119931 ) , iBode et al.l (11994 ) , iGarijo. Athanassoula. fc Gc 

(119971 ). These simulations find a runaway merging process occur ing near the cluster c enter 

forming a BCG. Recent analyses of the Millenium Simulation (jSpringel et al.l (120051 )) by 

Lucia et al.l (120061 ) indicate that, on the one hand, the most massive elliptical galaxies have 

the oldest and most metal rich stellar populations. On the other hand, massive ellipticals are 

predicted to be assembled later than their lower mass counterparts, and they have a larger 

effective number of progenitor systems (reaching up to A^^efr ~ 5). 
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In order to investigate the influence of galaxy mergers on the BCG mass dispersion sys- 
tematically and with statistical significance, we set up ensembles of Monte Carlo simulations 
of cluster evolution. A set of parameters that describe the initial conditions of the cluster 
is defined, then the cluster evolution is simulated and the evolved cluster is characterized. 
After that, the simulation is repeated from the beginning with the same initial parameter set, 
evolving to a different cluster due to the somewhat stochastic nature of galaxy encounters. 
The number of simulations of cluster evolution for a fixed set of initial conditions defines 
the ensemble size, which in turn determines the level of confidence with which parameters 
of the evolved cluster can be specified. Then one of the initial cluster parameters is varied 
leaving all other parameters fixed, and an ensemble of cluster simulations is run in the same 
way as described above. This parameter variation is repeated with a certain step size within 
an interval given by observations. The interval is supposed to cover all occurences in na- 
ture (e.g., the cluster richness) and/or uncertainties as the true value could not yet be fixed 
by observation s (e.g. , the fraction of cluster mass that is still bound to individual galaxies. 



Springel et al.l (120011 )). Having gone through all steps of the variation of this particular 
cluster parameter, the next parameter is varied leaving all others fixed. In this way, the 
full parameter space is covered, while all parameters are varied independently of each other. 
Possible correlations of initial cluster parameters (e.g., size and richness) can easily be ap- 
plied afterwards by restricting the parameter space to a subspace. The required computing 
power for such a large number of repeated simulations of cluster evolution over a Hubble 
time, covering a broad parameter space, is still beyond the capacities of present-day detailed 
N-body simulations. That is why a Monte Carlo code has been developed t o investigate th e 
small BCG mass dispersion, surpassing a purely analytical description (e.g., Merritt ( 19831 )) 
with the power of simulating statistically significant ensemble sizes. 

The structure of this paper is as follows: In the next section the simulation method is 
described. Section 3 details the parameter set and the mass dispersion. Section 4 presents the 
results of three simulation sets with exemplarily representative parameter sets. A discussion 
of the results is given in Section 5. Finally, Section 6 briefs the conclusions of this work. 



2. The model 

Numerical simulations require the masses of the galaxies. On the presumption of a given 
mass-to-light ratio T, Eq. ([T]) is directly translatable into the BCG mass dispersion. Though 
masses resulting from dynamical simulations are still dependent on T once compared with 
the observed luminosities, the relative variation does not depend on T. 

The following subsections briefiy describe the model assumptions that enter the simu- 



lation computations. For further details see iJagemann 



fll997t ). 



2.1. The cluster 



The cluster simulation starts with galaxies homogeneously distributed in a sphere of 
radius R and total mass M. R and M are kept constant in time in the course of development; 
in other words, we assume that the cluster is isolated and decoupled from the Hubble flow. 
At the beginning the number of galaxies in the cluster equals N = Nq. N diminishes in 
the course of time due to mergers and destruction of galaxies. The mass of the destroyed 
galaxies is counted to the intergalactic medium, so that the total mass of the cluster remains 
constant. 

The galaxy masses are initially drawn from the exponential distribution function dis- 
cussed in Sect. O The model starts from the assumption that the total mass of the galaxy 
cluster is proportional to the total mass of all galaxies: 



Mtot = AroMo-fc (2) 

Here A^o is the initial galaxy number and Mq is the mean galaxy mass (Sect. [3]). Then the 
constant k is the ratio of cluster mass to galaxy mass, i.e. the ratio of mass-to-light ratio of 
the cluster (in the following denoted by Tdustor) and mass-to-light ratio of the galaxies: 

7 ^cluster 



T 

galaxies 

The total mass of a (rich) cluster is of the order of IO^^Mq. Therefore the mass-to- 
light ratio of a cluster is within the range [200; 500] Tq, and so the discrepancy between the 
mass-to-light ratios of present clusters as a whole and their individual members is at least 
k = 5 and probably more. An interesting and crucial question concerning the interaction 
between galaxies of a cl uster is h ow m uch of the cluster dark matter is really gravitationally 



bound within galaxies (iMerrittI (Il985l )). i.e. the value of k. Two possible extreme cases are 
imaginable: either the dark matter hides within individual halos of single galaxies, or these 
halos are smeared over the whole intergalactic background due to tidal processes or similar 
effects. As virial velocities of the galaxies are determined by the total cluster mass, whereas 
individual galaxy masses enter decisively those equations that determine the interaction 
between two encountering galaxies, the trend resulting from the different extreme cases can 
be anticipated: 



-6- 



If dark matter exists in individual halos, the merger cross section is much larger than 
in a cluster containing significantly intergalactic matter, so the model develops more rapidly 
and more violently. This work considers these two cases. 



2.2. The galaxies 

In this paper individual galaxies are described by two parameters: their radius rj and 
their mass Mi (i — 1 . . .N). Mass and radius are changing in the course of cluster evolution; 

at the beginning, the masses of cluster galaxies are characterized by a distribution function 
with characteristic mass Mq and an associated distribution function in galaxy radii with 
characteristic radius Tq. 

The evaluation of galaxy encounters is done under the assumption that galaxies are 
represented by Plummer spheres: 

^ 3_ Ma- 

Here a is the Plummer radius, M is the total mass of the galaxy and G is the gravitational 
constant. Later, for the energy transfer in the tidal approximation, the root mean square 
(RMS) radius of the extended perturbed mass distribution is needed. So we further assume 
a cut-off radius trms ~ 4.3 a. In the following, we call this proportional factor /3 := rRMs/o- 
The galaxy radii n = trms are related to the drawn masses via the Faber- Jackson-law. 
Calculating interaction of cluster galaxies, the significant parameter is not the radius but 
rather the binding energy which computes for a Plummer model to: 

E=-^^. (5) 

64 a ^ ' 



2.3. Encounters 

Encounters of galaxies let the model evolve. In this work, encounters are described by 
the two parameters v and B. v is the relative velocity, and the maximum impact parameter 
B defines the maximum separation between two passing galaxies that is evaluated to be an 
encounter. Fixing the total cluster mass M^ot and characteristic radius R to be independent 
model parameters, the mean square of the spatial velocity is set by the virial theorem 
that is assumed to be apphcable to the cluster: 
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2 ^GMokN 
^3 = / ^ (6) 



This procedure is in contradistinction to the procedure of iKrivitsky fc Kontorovich 



(Il997l ). They find a kind of "explosive" merging without the additional constraint of virial 
equilibrium. 

The radius R of the sphere is assumed to be proportional to the gravitational radius 
with proportionality factor In addition, we assume the spatial virial velocity fs (that 
also determines the rate of encounters) equals a/ ("W"^) that has to be put into the formula 
of energy transfer. 

The impact parameter h of the actual encounter is chosen within the limits of [0; B] such 
that encounters are uniformly distributed upon the maximum encounter cross-section TtB"^. 

Using the impulse approximation, in this model B is defined the way that the error due 
to neglecting encounters with impact parameters greater than B is less than one per cent 
for a galaxy of standard size tq [tq corresponds to Mq via the Faber- Jackson-law). 

The mean time between two encounters is given by 

. _ 1 
^ ~ 2>vB^ N{N - 1) ^ ' 

and diminishes with decreasing galaxy number during the cluster development. The total 
development time is set to T = 10^°yrs. 



2.4. Development 

Gravitational interaction is the origin of cluster development. We compare the kinetic 
energy from orbital motion of two encountering galaxies with their internal energies. The 
binding energy of an individual galaxy increases at the expense of the orbital energy, i.e. 
at the expense of the total kinetic energy of their relative motion. As long as the total 
energy transfer AEi + AEj is less than the kinetic energy of the relative motion, the pair 
of galaxies is assumed to stay unbound. The energy input leads to an increase in binding 
energy of a galaxy which therefore expands. If energy input exceeds the binding energy of 
the concerning galaxy, this leads to its destruction. 

But if the total energy transfer is higher than their kinetic energy. 



AE, + AEj > Ekin, 



(8) 
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the encounter ends up in a bound system, and it is assumed that the galaxies merge. If the 
encounter leads to a merger event, both participating galaxies are replaced by a single one 
(the remnant), and its mass and radius has to be recalculated. The remnant mass is simply 
the sum of both galaxy masses, and the radius is computed from energy conservation: 



Er = Ei + Ej + E] 



kin' 



(9) 



It is possible that the remnant binding energy becomes Er > 0, if 



-E'kin > —{.Ei + Ej). 

In this case no remnant is created and both galaxies break up. 



(10) 



To calculate ene rgy exchange, the energy transfer is expressed using the impulse and 
tidal approximation (jSpitzerl (119581 )) combined with a smooth interpolation of the two lim- 
iting cases of 6 = and h>h{h = 5max(rh,i, Thj), is the median radius of the respective 
perturbed system). 

The energy transfer of a particle of mass Mj passing an extended particle system at 
distance h with velocity v according to impulse and tidal approximation is 



AG'^M'jMirf 
= TTki • 



(11) 



where rj is the RMS-radius of the perturbed system i ( ISpitzen (119581 )). 



In the case of central penetration 
Plummer models is given by 



0) the energy transfer between two identical 



(12) 



where (3 was introduced in Subsect. 12. 2[ The interpolation formula is then assumed to be 



Ai^^ma; 



l + ft&^ 



(13) 



(Fig. p. 



Proceeding this way we use two terms of energy transfer in high speed encounters, i.e. in 
encounters where relative velocities are higher than the galaxy internal velocity dispersion. 
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Fig. 1. — Energy transfer with interpolation formula Eq. [131 Abscissa: 6/6, ordinate: 
A£'/A£'inax; dashed line: impulse approximation, solid line: interpolation. 
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This approximation is not too crude (lAguilar fc Whitd (Il985l )) although galaxy relative 



velocities are about as high as their star velocities 

Other dynamical processes, especially dynamical friction between galaxies or galaxies 
and the intergalactic background, tidal interaction between galaxies and the cluster potential 
as well as mass-loss due to ram-pressure are neglected. 



3. Mass dispersion 



Trying to understand the BCG mass dispersion, it is necessary to know more about 
the stochastic variance, inherent in the distribution function. When drawing galaxies 
from a given frequency distribution, there is already a dispersion in the mass of a distinct, 
e.g. the brightest, galaxy depending on the kind of distribution and the number of drawn 
galaxies. An analytical approxim ation of th e glob al luminosity distribution of galaxies in 
today's clusters is Schechter's law (jSchechteii (119761 )): 



N{L) dL 



AT. (A 



with 



dL 
J Lit 



(14) 



a = -5/4 

= 1.0 ■ lO^'' h~'^ Lq in the visual band 
A^, G [20; 115]. 



a is given within an error of 20%, but there are also votes for an exponential distributio n 
in compact clusters (a = 0) differing from that of field galaxies (iGeller fc Peebles I (119761 )). 
In any case, solely the number of th e brighte s t (cD -) galaxies in the centres of galaxy clus- 
ters are not exactly accounted for ( Dressier (jl978l ). and references therein). This fact by 
itself is an indication of the answer to the central question: Is Schechter's law the result of 
development in kind of mer ger and destr uction processes or does it reflect the mass spec- 
trum in the early universe (ILauerl (119851 ))? Since the BCG is in an exceptional position 
with respect to Schechter's law and does not fit the global distribution, it is not probable 
that Schechter's distribution accrues from a different distribution (except that other than 
development processes would affect the BCG luminosity). Otherwise the distribution would 
include the BCG. 



With the assumption of a specific mass-to-light ratio of individual galaxies, galaxy 
masses can be drawn from a universal mass distribution. Because of better analytical han- 



dling and to avoid divergences galaxy masses of the model cluster are drawn from an expo- 
nential distribution. In this case, the initial probability distribution is 



Pbcg(M) dM = 1 - e «o 




AL\N-1 1 




M 



(15) 



shown in Fig. [2l 

In this case the expected value of the most massive initial galaxy can be written as an 
harmonic sum that depends on the mass Mq characterizing the distribution function and on 
the number of drawn galaxies Nq: 



With = 1000, e.g., the relative deviation computes to (Tbcg/ (^bcg) = 0.17. Though 
the expected value of the BCG masses increases with galaxy number, the relative deviation 
does not. Generally the fact is remarkable, that initially BCG masses depend on the number 
of galaxies, whereas in reality such a dependence is not found, and the BCG relative mass 
dispersion is distinctly smaller for all numbers of galaxies than the observed one. These 
differences may be due to development processes, although additional correlations among R 
and may cause selection effects in {N, i?)-space. 

Statistical statements about clusters of galaxies are derived from model assumptions 
depending on several parameters. Thus the resulting quantities to be analysed are only 
estimates of the mean sampled over infinite runs. The accuracy, given as relative la variance, 
is proportional to the square root of the number of simulation runs. Differently set up model 
clusters with the same set of initial parameters are here called an ensemble. So since all our 
simulation results are based on ensemble sizes of 100 runs, they are accurate to within 10%. 

Additionally, any dispersion in absolute magnitudes (or masses, respectively) in a model 
that depends on variable parameters will generally depend on these parameters, too. For 
example, BCG masses - as well as their dispersions - that are drawn from a distribution 
function depend on a characteristic mass Mq scaling this function and depend also on the 
total number of galaxies A^o drawn from this distribution. 




(16) 



k=l 



The BCG mass dispersion is then computed to 




(17) 
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Fig. 3. — BCG masses and their variance at the beginning of the simulation. Initial param- 
eters are: Mq = 4 • lO^^M©, ro = 15 kpc, vq = 800 km/s, R = 3 Mpc, T = 10 Gyr. 
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Fig. 4. — BCG masses and their variance at the end of the simulation. The initial parameters 
are the same as in Fig. [31 
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Figures [3] and m illustrate the dependence of the mean of BCG masses on the initial num- 
ber of galaxies. Likewise, standard deviations belonging to any value of N are shown, too. 
Figure [3] describes the initial state. Fig. IHthe evolved state. The curve in Fig. [3] is smooth, 
because the expected value of BCG mass and dispersion are theoretically computable. In 
Fig. [3] the stochastic 10% noise of the mean values is obvious. 

In order to calculate the total dispersion of BCG masses, let these masses be a function 
of a vector x in parameter space. Furthermore, assume that a distribution function pg{M) 
is given. The expected value of M then averages the ensemble and the parameter space to 



((M)) := {M{x))p{x)dx, (18) 

J X 

where {M{x)) means an ensemble average. The mass dispersion averages the ensemble and 
the parameter space to 



(AM)L = {{AMY) + {{M)') - ((M))^ (19) 

where the inner brackets always denote an ensemble average. Thus, the total dispersion is 
the expected value of the individual dispersion in x plus the scatter of the expected values 
of M in x: 



(AM)L = {{AMY) + {A{M)f. (20) 

This expression is used later on to calculate the mass dispersion of the simulations. At the 
beginning of the dynamical development, any scatter in BCG absolute magnitudes or masses 
is the sum of different kinds of model dispersion: 



1. The statistical variation at constant model parameters within an ensemble, whereby 
one gets a convolution of two distributions: the stochastic variance at a constant vector 
of model parameters and 

2. the variance caused by the distribution of clusters upon the model parameters. 



Hence the observed relative mass dispersion is the upper limit of the statistical error. 



In this context it is astonishing that there is definitely no correlation be tween Mrcq an d 
richness (meaning the fi nal galaxy number) in today's observed clusters (jSandagd (119721 ): 



Postman fc Lauerl (119951 )). But drawing galaxies from a distinct distribution (except from a 



5-function), there is always - within the bounds of resolution - such a dependence, so that 
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there inevitably has to be another process to explain the independence of Mbcg on N. The 
simulations try to answer the question whether mergers can account for such a behaviour. 
The observed scatter in absolute magnitudes is in any case an upper limit of the intrinsic 
BCG mass dispersion within any ensemble. 

The dependence of the BCG absolute magnitudes on the scaling mass Mq is more 
complicated. Mq is positively correlated with Mbcg with high statistical significance. In 
order to eliminate any photometric calib r ation error as well as a global shift of the luminosity 



function, iTrevese. Cirimele. fc Appodial (Il996l ) compare Mbcg ~ ^lo (Mio is the magnitude 



of the 10th brightest galaxy) with Mq — Miq. They find, that both differences are negatively 
correlated. 



4. Results 

Table [T] lists the input parameters of the simulation. 

Here the two parameters R and are allowed to vary independently (2) in the given 
interval while the pair (Mq; tq) remains constant. Then, the mean relative velocity computed 
using the virial theorem adjusts automatically. Afterwards a new pair (Mq; tq) is fixed from 
the variation interval, and A^ and R are varied again. Mq and vq are not correlated, k is set 
to 1 or 10 (the probable value of present clusters). 

The underlying idea is that there is a fixed value of Mq for all clusters in nature within 
narrow limits (corresponding to the fixed in the Schechter distribution), that the mass- 
to-light ratio of all clusters is initially the same and that all galaxies of equal initial mass 
are of the same size, independent of cluster membership. 

On the other hand, there are clusters of different size and variable richness in nature 
that enter the measurements of mass dispersion of the brightest cluster galaxies and are 
equally put together. 

At the end of the simulation the following ensemble mean data are derived: 

1. The mean relative velocity calculated from the virial theorem, 

2. the number of galaxies at the end of the simulation, 

3. the total number of encounters within the maximum impact parameter, 

4. the number of constructive merger processes. 
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5. the number of destructive merger processes, 

6. the number of destroyed galaxies without forming a remnant, 

7. the number of merging events that lead to the formation of the BCG, 

8. the BCG mass, 

9. the BCG mass ensemble dispersion. 

Except for the last two items giving the relative mass dispersion, the data serve as 
investigation of the model and are not connected to observations. 

4.1. The simulation with Mq = 5 ■ lO^^M©, ro = 20kpc and k = l 

In the following, results of this parameter set are discussed in detail. Concerning other 
values of Mq, tq and k only occuring changes in contrast to this set are mentioned. 

The choice of Mq = 5- IO^^Mq means a moderate mass-to-light ratio of about T galaxies = 
10 of today's observed clusters, i.e. the absence of dark coronae, and k = 1 means the initial 
absence of intergalactic matter. 

The number of encounters of all galaxies is shown in Fig. O This total number of 
encounters increases within the bounds of parameter variation to 10^ and corresponds es- 
sentially to Eq. [71 substituting the total development time T for At. Therefore the number 
of encounters is proportional to N^^"^ and to i?"''/^ and changes in the course of evolution 
only at small radii R , where the galaxy number decreases due to merger and destruction 
processes (cf. Fig. ITUi). 

Figure [6] shows the mean BCG masses at the end of the simulation. They are located 
between about 30 and 70 ■ lO^^M© and show a clear peak at small radii and small initial 
galaxy numbers. 

In order to exclude the influence of the initial mass distribution in this histogram, the 
pure mass increase is shown in Fig. [71 Generally the BCG is growing in the course of 
development. The increase amounts to 40 -lO^^M© at the peak which stands out even more 
clearly than in Fig. [6] 

The main result of the simulation, i.e. the mean values of the BCG relative mass 
dispersion, is shown in Fig. [8l The dispersion is nearly constant at about thirty per cent and 
increases at the small clusters to a maximum of 57 %. 
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Table 1: Simulation parameters 



Parameter 


Value (s) 


Remark 


Ensemble size 


n = 100 


constant 


Total time 


T = 10 Gyr 


constant 


Char, mass 


Mo e {5; 50} • lO^M© 


variable (1) 


Char, radius 


ro e {20; 100}kpc 


variable (2) 


T ratio 


k e {1;10} 


variable (1) 


Cluster radius 


R e [0.25; 3.00] a 0.25 Mpc 


variable (2) 


Richness 


N e [100; 1000] a 100 


variable (2) 



If 




Fig. 5.- 



Logarithm of the total number p of encounters. 
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Fig. 8. — Mean values of BCG mass dispersion after the simulation, d — (ctbcg) / (Mbcg) ■ 
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Fig. 9. — Mean values of the increase of BCG mass dispersion. d 

(cBCG.afterwards)/ (-^BCG .afterwards ) — (cBCG,previously)/ (-^BCG ,previously/ ■ 



- do ^ 
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Fig. 10. — Mean number of mergers leading to the BCG formation. 
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In order to investigate the increase in mass dispersion as in Fig. [7] and to exclude the 
influence of the initial distribution, the relative mass distribution after the simulation is 
subtracted by the relative dispersion before the simulation in Fig. O As in Fig. [TJ the peak 
stands out more clearly against the stochastic noise. 

To answer the question how many mergers are necessary to form the BCG this parameter 
/ is shown in Fig [TOl The maximum value lies at about six mergers and decreases steeply 
from there to a value of about one. 

Not only merger events but also destructions of galaxies are part of cluster development. 
Figure [11] shows how many galaxies are left at the end of the simulation. Here a clear trend 
is recognizable: the percentage of destroyed or merged galaxies increases considerably with 
decreasing cluster radius. In the case of rapid evolution the number of galaxies may decrease 
down to twenty per cent of the original value. 

The difference between Fig. [11] and Fig. [13] provides roughly for the number of galaxies 
destroyed in encounters (Fig. [T2|). It turns out that the biggest fraction of galaxies is de- 
stroyed when the energy transfer is less than the sum of binding energies of the 2 galaxies. 
The percentage of disrupted galaxies is always less than 65 per cent and decreases with in- 
creasing radius. At small radii the fraction of destroyed galaxies is roughly independent of 
the initial galaxy number, but increases with galaxy number at increasing radii. 

Finally, to address the question about the fate of the remaining galaxies, the number of 
all the merger events of all galaxies is divided by the final galaxy number in Fig. [13] The 
parameter t in this diagram denotes the percentage of those galaxies that are involved in a 
merger event. In the same way t is a measure of the frequency describing how often a galaxy 
has been involved in a merger event averaged over all cluster galaxies. In the peak vicinity 
every galaxy has experienced one merger on the average, t sharply decreases with increasing 
and R. 



4.2. The simulation with Mq = 5 ■ IO^^Mq, tq = 20kpc and k = 10 

The choice of Mq and tq corresponds to the choice in Sect. 14. H A; = 10 corresponds (at 
Mq = 5 ■ IO^^Mq) roughly to the maximum value of the present fraction of intergalactic dark 
matter in galaxy clusters. 

The mass distribution of single galaxies corresponds exactly to the one of Sect. 14.11 but 
because of the addition of intergalactic matter the virial velocity is increased by the factor 
vTO according to Eq. [H] This is why the number of encounters is about this factor higher 
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Fig. 11. — Mean number of destroyed/merged galaxies per initial galaxy number after de- 
velopment. 
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Fig. 12. — Percentage of destroyed galaxies. 
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Fig. 13. — Mean number of merger events per final galaxy number. 
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in comparison with Fig. [5l The high spatial velocity-dispersions a little above the diagonal 
of more than 2400 km/s are not observed in reality; therefore this part of the simulation is 
excluded in the following discussion. 

The final galaxy numbers drop to values of significantly less than hundred per cent 
at small radii. The qualitative behaviour follows the one of Fig. [TTl but the maximum 
is, however, clearlier marked at small radii and galaxy numbers, and about half as many 
galaxies less have been destroyed at a fixed pair {R,N) compared with Sect. 14.11 The 
dominant development process here is not merging but destruction of encountering galaxies. 
The mean mass distribution of BCGs after the simulation has not changed much: The 
initial mass distribution can easily be recognized. Deviations in comparison to the initial 
distribution arise in about half of the cases, these deviations are about ±1 ■ lO^^M© for each 
half. 

Only little signs of evolution are present in the BCG mass dispersion after the simulation: 
The initial mass scatter is still present. The relative dispersion does not exceed a value of 25 
per cent. The trend of the initial distribution having smaller spreads with increasing galaxy 
numbers still exists. Just as for the mean mass values there is an irregular change in the 
scatter of 1-2% compared with the initial scatter. 

4.3. The simulation with Mq = 50 • IO^^Mq, tq = lOOkpc and k = l 

The total mass of the galaxy cluster is the same here as in Sect. 14. 2[ the galaxies however 
are ten times more massive, k = 1 means the initial absence of intergalactic matter as in 
Sect. 14. 1[ The dark matter is now hidden in individual halos instead of spread throughout 
the intergalactic medium. Therefore more extended galaxies with a characteristic radius of 
To = 100 kpc are chosen at this parameter set. 

Figure [13] shows the BCG mass mean values after the development. In comparison with 
Fig. E] the shift of the maximum toward smaller galaxy numbers but significantly greater 
cluster radii is striking. 

Figure [15] makes it clear that the relative mass dispersion remains always below forty 
per cent and even decreases significantly going to greater galaxy numbers, particularly at 
this set of parameters. 

Fig. [16] shows that the number of mergers leading to the formation of the BCG is non- 
zero for a wider range of parameters N and R. The maximum mean number of mergers is 
four as in Sect. 14. 1[ 
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Fig. 14 



, — Same as Fig. El but now for model 14.31 
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Fig. 15. — Same as Fig. [HI now for model H73l 
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At smaller galaxy numbers the percentage of those galaxies that have been involved in 
merger processes increases up to sixty per cent (Fig. [T8|) . cf. also Fig. [131 From a galaxy 
number of about four hundred onwards there are practically no merger events to be noted. 



5. Discussion 



There are two processes with opposite effects which rule the evolution of the model: 
On the one hand the number of encounters decreases with increasing cluster radius at fixed 
galaxy number (cf. Fig. [5]). On the other hand the virial velocity increases with increasing 
total cluster mass by increasing richness at fixed cluster radius. Although the number of 
encounters rises with velocity, the equations of energy transfer (Eqs. [11] and [12]) depend 
on and the kinetic energy term of the relative motion depends on . So high energy 
transfers are only effective for mergers at low velocities. Taking both effects together, only 
compact groups evolve noticeably in contrast to clusters. The number of encounters g oi 
single galaxy for the total development time can be estimated to 



if richness does not change dramatically. The merger probability can be estimated in the 
same way: let hm be the impact parameter so that the merger criterion is fulfilled within its 
enclosed area: 



= (22) 

Then the fraction of merger events (caused by the stochastic choice of impact parameters) 
is 6^/5^. This fraction multiplied by the number of encounters then reveals the number of 
merger events V of one galaxy (e.g. the BCG): 



V = %r^r^T^G^^-s^— (23) 

So the number t of mergers depends strongly on the cluster radius and decreases with 
increasing R (cf. Fig. [TT]l . On the other hand t increases proportional to the square root of the 
initial number of galaxies (assumed that it remains nearly constant, e.g. Fig. [lT]at R = 0.75 
Mpc and higher). The boost of constructive (remnant forming) or destructive encounters 
can thus be explained, but not the merger boost at low galaxy numbers. Furthermore it has 
to be taken into account whether the remnant remains bound. That leads to 
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Fig. 17.— 



Same as FigJTT| now for model 
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Fig. 18. — Same as FigHSl now for model 1131 
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^ > 1. (24) 



8ro/2 N 

The remnant has a chance to survive only at high cluster radii and small galaxy numbers. 
Together with the merger criterion this leads to the occurring maximum development at small 
cluster radii and small numbers of galaxies (Figs. [71 [TUl [T51) : Small radii are necessary to 
form remnants, and their binding energies are negative only at small galaxy numbers. 

These estimates expressed in simple formulae (Eqs. [2T1I231) are possible because our 
simulation show that there is no agravating change of the initial mass distribution. 

The simulation computations of Sect. 14.31 show a merger behaviour of all galaxies and 
especially of those forming the BCG comparable to Sect. 14. 1[ A higher galaxy mass has a 
positive effect on the number of mergers because the mean energy transfer rises with the 
third power of the characteristic mass, whereas the kinetic energy of the relative motion still 
rises linearly. On the other side a higher cluster mass in Sect. 14.31 in contrast to Sect. 14.11 
causes the virial velocity to increase, suppressing energy transfer such as in Sect. 14. 2[ 

Interpreting simulation 14.21 to represent cluster states further evolved in time compared 
with sirn u lation 14.11 as the mass of the intergalactic medium increases due to tidal stripping 



(IMerrittl (119841 )). a picture of cluster evolution is formed: In the early phases of cluster 



evolution merger processes dominate the d evelopment of m ass spectrum and BCG, then 



paling into insignificance as time goes by (cf. lAlbingerl (119951 )). Comparing both simulations 



it should be noted that the cluster mass is the same in simulation 14.21 at = 100 and 
simulation 14. II at A^ = 1000. An increase of the mass fraction of intergalactic matter leads to 
a drastic decrease in the merger rate at equal galaxy mass distributions, because the mean 
kinetic energy controlled by the virial theorem (Eq. [6]) outgrows the expected energy transfers 
in Eq. [131 so that mergers become improbable. The amount of galaxy binding energy is not 
influenced by a variation of fc, and so a (less steep) decrease in galaxy destructions is expected 
on the basis of Eq. [T3] as shown by the simulations. 

A closer consideration of Fig. [7] together with Fig. [TOl dividing the mean increase of 
BCG mass by the mean number of mergers leading to the BCG formation, shows that the 
brightest galaxy must have been built by mergers of galaxies of uniform mass Mq independent 
of general con ditions of the cluster. This is in accordance with the observation results of 



Hoessell (119801 ). 



Figure [13] gives information about the number of constructive mergers per final galaxy 
number, i.e. the fraction of galaxies that merged with other galaxies during the simulation. 
This percentage is equivalent to the averaged number of mergers a galaxy has undergone. 
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This fraction rises going from high to smal l clust er radii and amounts to between some 



and about hundred per cent. Toomre (IShul (119821 ) ) speculated that perhaps all elliptical 



cluster galaxies have been built by mergers. From the fraction of galaxies that are presently 
involved in gravitational interactions and the probable life time of this passing phenomena 
he calculated that merging of spiral galaxies happens for a fraction of about ten per cent of 
all cluster galaxies, in good agreement with the results of this work. The fact that especially 
rich compact clusters which are assumed to be dynamical older systems than spiral rich ones 
have a signif i cantly larger fraction of elliptical galaxies (up to 35 per cent in the Coma cluster 



(cf. lOemlerl (Il974j ))) may have multiple reasons. First a comparison of the simulations in 
Sects. 14. II and [4.21 shows that mergers grind to a halt at rising fraction of intergalactic matter, 
as the fraction of prospective interacting galaxies decreases at constant total cluster mass 
due to mergers and destructions. Hence the merger rate is higher in a cluster of earlier 
state of development (corresponding to irregular spiral rich clusters) than in the final state 
of compact clusters. So the fraction of elliptical galaxies in Toomre's estimation has to be 
corrected upwards. Secondly elliptical galaxies may be a result of disturbances caused by 
nearby passages without merging that would lead to its rise to the presently observed value. 

The percentage of galaxies destroyed in encounters (Fig. [T^ yields an estimate of the 
dark matter fraction in clusters provided by destroyed galaxies. This fraction amounts up 
to 60 per cent corresponding to a ratio of cluster mass and total galaxy mass of k = 2.5. 

Further simulations with the presented model were carried out using other sets of pa- 
rameters without leading to different results than discussed in this work. Particularly an 
increase in galaxy numbers beyond a value of 1000 does not change anything, because the 
evolution of any model with scales comparable to clusters takes an imperceptible course at 
high galaxy numbers. 



6. Conclusions 

There is a common result from all simulations of Sect. HI the cluster evolution is tightly 
related to the evolution of the brightest galaxy of this cluster. Only at small galaxy num- 
bers and simultaneously at small cluster radii the dynamical development due to the model 
described in Sect J2] plays a noticeable role (cf. e.g. Fig. [10]). At the same time increasing 
merger action is a concomitant of increasing destruction of galaxies. So, in the parameter 
range of rapid evolution, the number of galaxies depleted such that the evolved small cluster 
looks like a compact group rather than an Abell cluster. Nevertheless merging is always the 
dominant process (at least for the most massive galaxies). For clusters of presently observed 
richness and their most luminous member (being the comparative value to our simulation) 
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the dynamical development of the mass spectrum and of correlated quantities as the BCG 
mass causes only little changes with respect to the initial conditions. 

Considering Fig. E] and Fig. [S] the steadiness of BCG masses and their ensemble disper- 
sion attract attention (irrespective of compact groups). This approximate independence of 
masses on cluster radius is observed in reality, and there are two causes that balance each 
other: At the beginning of the simulation the BCG masses are small and rise with increasing 
richness. In the course of evolution mergers take place preferentially at low galaxy numbers 
resulting in a relatively constant BCG mass level after the development in dependence of the 
parameters R and N. Tak i ng in to account an additional correlation between cluster radius 



and richness (e.g. lOemlerl (Il974j )) the BCG mass steadiness would turn out almost clearer, 
just as their dispersion would be narrower than stated in this work because of the restriction 
in parameter space. 

Why does the BCG mass scatter remain significantly narrow even at the end of cluster 
development? First of all the number of mergers is relatively small for the mentioned reasons 
(between one and six merger events depending on the parameter set), and secondly the 
identity of a galaxy to be BCG is fixed only when the cluster development is complete. 
So, not the mass dispersion of a distinct galaxy fixed before the starting development is 
investigated, but the mass dispersion of a galaxy that is a posteriori identified to be the 
most massive one. Therefore the BCG mass dispersion turns out to be smaller than in the 
case of any individual galaxy chosen at the beginning of the simulation and following its 
evolution. 

I thank Prof. Dr. R. Wielen and Prof. Dr. B. Fuchs for stimulating discussions and 
Dr. S. Frink for help in the preperation of this work for publication. 
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